Nonlinear dynamics of a microswimmer in Poiseuille flow 
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We study the three-dimensional dynamics of a spherical microswimmer in cylindrical Poiseuille 
flow which can be mapped onto a Hamiltonian system. Swinging and tumbling trajectories are iden- 
tified. In 2D they are equivalent to oscillating and circling solutions of a mathematical pendulum. 
Hydrodynamic interactions between the swimmer and confining channel walls lead to dissipative 
dynamics and result in stable trajectories, different for pullers and pushers. We demonstrate this 
behavior in the dipole approximation of the swimmer and with simulations using the method of 
multi-particle collision dynamics. 

PACS numbers: 47.63.Gd, 47.63.mf, 47.61.-k 



Microswimmers often have to respond to fluid flow and 
confining boundaries, like sperm cells in the Fallopian 
tubes [lj| or pathogens in blood vessels Q • Artificial mi- 
croswimmers constructed with the vision to act as drug- 
deliverers in the human body Q would have to swim in 
narrow channels like arteries. Two properties influence 
the swimming in microchannels. On the one hand, vor- 
tices in flow reorient the swimming direction of microor- 
ganisms. In simple shear flow, for example, microswim- 
mers tumble due to a constant flow vorticity Q ■ Vortices 
in Poiseuille flow in combination with bottom-heaviness 
due to gravitation lead to stable orientations of swim- 
ming algae cells On the other hand, microorganisms 
swimming near surfaces are trapped by hydrodynamic 
interactions Q and ultimately escape with the help of 
rotational diffusion (Io| . Finally, bacteria in Poiseuille 
flow show a net-upstream flux at the walls due to the 
interplay of confinement and flow vorticity 0, H|. All 
these examples show there is genuine interest in under- 
standing generic features of microorganisms and artificial 
swimmers in Poiseuille flow. 

In this letter we demonstrate that the dynamics of a 
simple spherical microswimmer in a cylindrical Poiseuille 
flow can be mapped onto a conservative dynamical sys- 
tem with the Hamiltonian as a constant of motion. In 
analogy to the oscillating and circling solutions of a math- 
ematical pendulum, we discuss in detail the swinging and 
tumbling motion of the microswimmer in 2D and gener- 
alize them to three dimensions. Hydrodynamic interac- 
tions with the channel wall treated in the dipole approx- 
imation introduce dissipation and the microswimmer as- 
sumes specific stable swimming trajectories depending on 
its type as puller or pusher. 

We first introduce the geometry We consider a point- 
like microswimmer that moves with a constant intrinsic 
swimming speed vq in a cylindrical microchanncl where 
a Poiseuille flow is imposed. Using a cylindrical coordi- 
nate system (p,tp,z) with the coordinate basis (p, (p, z), 
the flow is given by v/ = — p 2 /i?Q h )z, where Vf is 
the maximum flow speed in the center of the channel 
[Fig. (Ha)]. In the absence of noise the equations of mo- 
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Figure 1. Swimmer in Poiseuille flow, (a) Flow profile 
v/(p), cylindrical coordinate system with coordinate basis 
(p, Cp, z) and orientation angle ^ for the projected orienta- 
tion into the p-z-plane. When e v — sin = 0, the motion is 
two-dimensional. Note that the sign of vorticity ft f changes 
when crossing the centerline. (b) Cross section of the mi- 
crochannel. The orientation in (^-direction defines the angle 
9. 



tion for the swimmer position r and orientation e are 
given by 



dt 



r = v e + v/, 



d 1 

— e = —Of x e (1) 

dt 2 1 V ' 



where $1/ = Vxvj = Vfp/R,Q b <p is the flow vorticity. 
The swimmer orientation e = e p p + e v ip + e z z has the 
components 

e p — — cos O sin *, e v = sin O, e z = — cos 9 cos 'J, (2) 

where \& € {— tt, it} is the angle in the p-z-plane [Fig.[TJa)] 
and 8 £ {— 7r/2, 7r/2} measures the orientation in the az- 
imuthal ^-direction [Fig. H{b)]. We note that 1*1 <7r/2 
means upstream and |*| > 7r/2 downstream orientation, 
respectively In the following we use rescaled units, 
p/Rch —> P £ {0,1}, z/Rch —> z and t/t — s- t with 
to = Rch/vo- We also introduce the dimcnsionlcss flow 
speed Vf =Vf/vo, which is the only essential parameter 
in our problem. 

First, we discuss 2D solutions of Eqs. ([T]) since they 
already capture many aspects of the swimmer dynam- 
ics. When O = 0, the trajectories of the swimmer are 
restricted to two dimensions, for example, to the x-z- 
plane, xe{ — 1, 1}. Due to the translational symmetry in 
z-direction, only the equations for x and * are coupled, 



2 



and Eqs. (JlJ give x 
results in 



sin*, * = VfX. Eliminating x 



* = vt sin *, 



(3) 



which is the equation of motion of the mathematical pen- 
dulum. Since in this analogy x plays the role of velocity, 
we can immediately write down the 2D Hamiltonian 



1_ 9 

iJ 2D = -VfX + 1 - COS* 

2 



(4) 



as a conserved quantity. Fig.[2]shows the a;-* phase space 
and typical trajectories z(x) for several flow strengths 
Vf. In analogy to the pendulum two swimming states 
exist. The flow vorticity rotates the upstream oriented 
microswimmer always towards the center. Hence, the 
swimmer performs a swinging motion around the cen- 
terline of the channel for i? 2 D < 2 which corresponds to 
the oscillating solution of the pendulum [e. g. blue tra- 
jectory of Fig. [2ja)]. For small amplitudes (* <C 1) the 
swinging frequency is ujq = ^/vj- When the upstream ori- 
ented swimmer moves exactly in the center of the chan- 
nel (stable fixed point), the Hamiltonian is zero. Down- 
stream swimming along the centerline (* = n) is an un- 
stable fixed point. After a slight disturbance of x = 0, 
vorticity rotates the swimmer away from the centerline. 
The swimmer performs a tumbling motion (i?2D > 2) 
which corresponds to the circling solution of the pen- 
dulum [green trajectory of Fig. [He)]. At H 2 d = 2, the 
separatrix x 2 = 2(cos* + l)/?;/ divides the swinging and 
tumbling region in phase space [red curves in the phase 
portraits of Fig. [2]. Since the Poiseuille flow is bounded 
by the channel walls, tumbling motion only occurs for 
Vf > 4 [Fig. [H[c)] . Meaning sufficiently strong vorticity to 
prevent the swimmer from crossing the centerline. 

If we only consider steric interactions of the swimmer 
with the channel wall, the swimmer crashes into the wall 
at \x\ = 1 for iJ 2 D >Vf/2, reorients due to the flow vortic- 
ity towards the upstream orientation, and leaves the wall 
at * = with H}^ x = Vf /2. The swimmer then performs 
a swinging motion between the walls with maximum am- 
plitude \x\ = 1 for Vf < 4 [green trajectory in Fig. Ha)]. 
So for vf < 4 the swimmer always enters a swinging mo- 
tion oriented upstream, at the latest after contact with 
the wall, whereas it tumbles close to the wall for Vf >4. 

To determine the full 2D trajectory in the microchan- 
nel, we solve the dynamic equation for z(t), 



z = Vf[l - x(t.) 2 } - cos*(i). 



(5) 



A careful analysis reveals the following. The swimmer al- 
ways moves upstream (i < 0), when Vf < 1—H2D, as shown 
in Fig.[2ja). When the flow is strong (vf > l+2H 2 u), the 
swimmer always drifts downstream (z>0), while swing- 
ing or tumbling [Fig. [2Jc)]. In between, mixed up- and 
downstream segments within one trajectory [Fig. HJb)] 
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Figure 2. Phase spaces a;-* (left) and typical trajectories 
z(x) (right) for several flow strengths vj. All trajectories start 
at z = 0. (a) upstream motion, (b) intermediate motion and 
(c) downstream motion. Note the various scales for the z-axis. 
The arrows indicate the orientation vector e of the swimmer. 



exist but a net upstream motion only occurs for Vf < 
I + i? 2D /2 [blue line in Fig.^b)]. 

For a non-zero azimuthal component, e v = sin ^ 
0, the swimmer trajectory is three-dimensional. Using 
Eqs. ((TJ) and Eqs. @, we obtain three coupled equations 
for *, © and p, 



p = — cos © sin * 
* = Vfp — sin tan© cos */p 
6 = sin6sin*/p. 



(0) 



Due to translational symmetry in z-direction and rota- 
tional symmetry about the channel axis, Eqs. ^ do not 
depend on z and <p. We are able to identify two constants 
of motion, 



L z = psin© 
H=\v f p 2 - 



I — cos * cos 9, 



(7) 



where L z is proportional to the angular momentum of 
the swimmer in z-direction. Due to this constant the 
sign of along a swimmer trajectory does not change. 
Eliminating 6 from Eqs. (JSJ), reduces the equations of 
motions to §§■ = — p, ^f" = v i / - So, H again plays the role 
of a Hamiltonian for the conjugate variables * and p. 

The intersection of the two constants of motion 
gives the orbit of the swimmer in p-*-© phase space 
[Fig. Eta, b)]. The stable fixed points of Eqs. ((6]) lie on the 
fixed-point line (p* = ^sin©* tan ©*/"/, ** = 0), drawn 
in Fig. Eta)- Swimming at a fixed point corresponds to 
a helical trajectory [Fig. [3Jc)]. The swimmer moves up- 
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Figure 3. p-SP-Q phase space. The intersection between L z — 
const, (orange) and H = const, (green) gives the phase space 
trajectory, (a) helical-like swinging motion (blue intersection 
curve) for L z = 0.2, H — 1. Black curve: fixed-point line 
corresponds to helical trajectories, (b) helical-like tumbling 
motion (blue intersection curve) for L z = 0.2, H — 3. (c-e) 
sketch of trajectories in the channel for helical motion (c), 
helical-like motion (d) and tumbling motion (e). 



stream for Vf < 1/cosO*, either on a left-handed he- 
lix (6 > 0) or a right-handed helix (9 < 0). Closed 
orbits around the fixed-point line correspond to swing- 
ing motion around a helical path [Fig. E[d)] and open 
orbits [Fig. E^b)] arc complicated tumbling trajectories 
[Fig. He)]. 

Now we consider hydro dynamic interactions of the mi- 
croswimmer with the bounding channel wall. The flow 
field for neutrally buoyant swimmers in a bulk fluid is 
in leading order a force dipolc, v(r) = [3(f -e) 2 — l]f 
where f = r/r and r\ is the viscosity of the fluid. For 
positive dipole strength, p > 0, the propelling apparatus 
of the swimmer is typically at the back (pusher), and for 
p < in the front (puller). Ref. [i[ treated the swim- 
mer close to a plane wall. The authors showed that hy- 
drodynamic interactions between the swimmer and the 
wall lead to re-orientation due to the wall-induced vor- 
ticity ttw and to attraction/repulsion. For example, a 
pusher swimming parallel to the wall is attracted and re- 
mains trapped at the wall, while a puller is repelled from 
the wall. We note that close to a wall the force-dipolc 
approximation for Jljy is no longer valid. Instead, ge- 
ometric details of the swimmer and thereby near fields 
become important for flw- Combining flow and wall ef- 
fects, the total angular velocity of a swimmer near a wall 
is fl(x 1 ^) = rtw(x,'^)+^if(x). Recent work showed that 
hydrodynamic re-orientation is almost negligible for bac- 
teria swimming near walls 0, H3] but rotates the swim- 
mer at the wall. For a sufficiently weak Poiscuille flow 
field, stable orientations ^>w for swimming at the wall 



exist when fl w (Rch,^w) + ^f(Rcb) =0 11]. Thermal 
fluctuations will, however, reorient the swimmer so that 
it leaves the wall [Toj | . 

In narrow channels a microswimmcr experiences hy- 



drodynamic interactions with the wall all the time. To 
capture the basic idea, we concentrate on 2D-trajectories 
and consider instead of the cylindrical channel wall, two 
parallel plates located at x = l and x = — l. Wc calculate 
the wall-induced translational and angular velocities us- 
ing the force-dipole approximation of Ref. [§[ and obtain 
the equations of motion, 



x = — sin 'F 
'F = VfX 



3p(3sin 2 <F - 1) / 1 



64tt 

3p sin \I> cos 'F 
64^ 



1 



(1-x) 2 (1 + x) 2 
1 1 



(1 - x)3 (1 + xf 



(8) 



where p = p/(rjvo) is the reduced dipole strength. Fig. [4] 
shows typical phase space plots generated from Eqs. dHJ 
for a puller (a) and pusher (b). The swinging motion 
of an upstream oriented puller becomes damped and an 
attractive fixed point in the center exists. Repelled by 
both walls, the puller swims upstream along the center- 
line. On the other hand, a puller tumbling near the wall 
is attracted, on average, by the wall and stays near to it. 
This is indicated by the green stable trajectory. All tra- 
jectories outside the unstable red limit cycle or separatrix 
converge to it. However, due to thermal fluctuations the 
puller may cross the separatrix. A pusher behaves differ- 
ently, it is attracted by the wall when oriented upstream 
in the center of the channel, but it is pushed away from 
the wall when tumbling near the wall. So, all trajectories 
converge towards a swinging motion about the centerline, 
characterized by a stable limit cycle in the z-'F-planc. 

To test our findings wc simulate the motion of a spher- 
ical microswimmcr in Poiscuille flow using the method of 
multi-particle collision dynamics (MPCD) 12J. It solves 
the Navicr-Stokes equations on a coarse-grained level and 
calculates the flow field around the swimmer in the cylin- 
drical microchannel taking into account both hydrody- 
namic interactions and thermal noise. In every simula- 
tion step randomly distributed point-particles of mass m 
at temperature fc^T first move ballistically for a time 
At and then they are sorted into cubic cells of length a. 
They interact with all other particles in the cell with a 
specific collision rule such that momentum is conserved 
locally. The density of the fluid is mN c / a 3 where iV c is 
the average number of particles per cell. Depending on 
the parameters and the specific collision rule, the viscos- 
ity r\ of the fluid can be calculated 12 , 3] ■ 

As a model microswimmer we use a spherical squirmer 
of radius R$ 13, [l^. It propels itself by a static, 



axisymmctric and tangential velocity field on its sur- 
face, v s (r s ,e) = (B 1 + (e-r s )B 2 ) [(e-f s )f s -e] where f s is 
the radial unit vector pointing from the center of the 
squirmer to the surface. The first mode B\ determines 
the swimming speed vq = 2Bi/?> and the second mode the 
strength of the force dipole p/rj = — AttB2R%, so B 2 > 
models a puller and B 2 < a pusher. The squirmer has 
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(b) pusher(p>0) 




Figure 4. Phase space trajectories for a microswimmer in a 
narrow channel for Vf = 10 for a puller (left) and a pusher 
(right). Green and red indicate, respectively, stable and un- 
stable trajectories, (a) and (b) are obtained from Eqs. JS} 
and (c) and (d) from MPCD simulations where we used the 
parameters B x = 0.045, /3 = B 2 /Bi = ±5, R s = 6, R C h = 18, 
N c = 30, At = 0.02 setting a = m = k B T = 1 and the colli- 
sion rule MPC-AT-a [l8j |. Each trajectory was obtained by 
averaging over 10 individual runs. 



already been used to investigate hydrodynamic interac- 
tions between several swimmers [lj, EEf and between a 
swimmer and a wall [16| and has been realized in exper- 
iments quite recently [17| . To implement the squirmer 
in the MPCD fluid, we follow previous work [TBI, [HI- 
Figs. Stc)-(d) show simulated trajectories for several ini- 
tial conditions for a puller (c) and a pusher (d). Our 
parameters are listed in the caption of Fig. 2J Although 
near field effects and the large extent of the swimmer play 
a role in the simulated dynamics, the qualitative behav- 
ior arising from hydrodynamic interactions between the 
swimmer and the cylindrical channel wall follows the an- 
alytical model. 

In order to learn about typical relaxation times in 
which swimmers approach their stable trajectories, we 
linearize Eqs. ((5J around the fixed point in the cen- 
ter. We obtain a harmonic oscillator equation for ^ 
with a friction coefficient 7 linear in the dipole strength, 
7 = — 3p/(647rr?i2p h ). For the squirmer, the estimated 
relaxation time for swinging motion becomes 7 _1 = 
to • 32^/(9/3), where (3 = B 2 /B u R = R C h/Rs, and 
to = Rch/vo w 1 s is a characteristic time scale for narrow 
microchannels. Typical values for (3 range from —1 to +1 
for existing microswimmers f2p| . So experiments should 
be able to observe that microswimmers approach their 



stable trajectories within seconds in sufficiently narrow 
channels. Similar estimates apply to the E. coli bac- 
terium where pm0.8pNu,m was measured recently [loT ]. 
Although hydrodynamic interactions between a small mi- 
croswimmer and a single wall may play no significant role, 
they become important in channels when the channel di- 
ameter is only few times the size of the swimmer. 

In conclusion, through a formal mapping onto a Hamil- 
tonian dynamical system we have shown that spheri- 
cal microswimmers perform either an upstream oriented 
swinging or a tumbling motion when moving in Poiseuille 
flow. Hydrodynamic interactions of the swimmer with 
the wall stabilizes the upstream orientation of pullers in 
the center of the channel whereas a pusher performs sta- 
ble oscillations around the centerlinc with a specific am- 
plitude. 

Spherical artificial swimmers with different locomotion 
mechanisms have been constructed and studied recently 
211 ]. Investigating them in microfluidic channels under 
Poiseuille flow, the generic features presented in this ar- 
ticle should be accessible in experiments. 

We thank R. Goldstein, I. Pagonabarraga, T. Pfohl, 
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Deutsche Forschungsgcmcinschaft for financial support 
through the research training group GRK1558. 
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